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SUMMARY 


^  j/A  theory  of  filamentary  structures  consisting  of  monotropic 
membranes  is  presented.  Applications  to  isotensoid  pressure  ves¬ 
sels  with  rotational  symmetry  demonstrate  the  use  of  the  theory. 
Particular  attention  is  given  to  applications  of  filamentary  design  of 
variable-geometry  expandable  structures. 


An  extension  of  the  theory  to  pressurized-filamentary  structures 
subject  to  centrifugal  loading  and  the  special  case  of  meridional^ 
winding  patterns  are  presented  in  two  Appendices. 


Physical  interpretation  of  the  resulting  shapes  and  winding  pat 
terns  leads  to  a  discussion  of  the  morphology  of  filament-wound 
pressure  vessels 


v — v 


Experimental  data,  obtained  from  filament-wound  toroidal  pres¬ 
sure  vessels,  confirm  the  validity  of  the  theory  and  demonstrate 
application  of  the  analytical  design  technique  for  filamentary  struc¬ 
tures.  „ 
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INTRODUCTION 


A  review  of  recent|  developments  in  the  field  of  advanced  struc¬ 
tures  for  aeronautical,  missile-borne,  and  space  applications 
shows  an  increasing  interest  in  exploiting  the  mirfaiduftfrfrr  phy sical- 
mechanical  properties  of  fibered  materials J  A  large  volume  of 
published  information  (See,  for  instance,  references  1,  2,  3,  and  4) 
is  available,  particularly  in  the  field  of  filament  wound  structures. 
^The  attraction  of  filamentary  structures  may  be  found  in  three  areasj^ 

/Td  Structures  of  / weight  ratios 

are  ^ate^i2i^£^possiblej(and  have,  in  certain  instances,  been 
demonstrated)  by  use  of  filamentary  materials  as  compared  to 
similar  structures  made  from  conventional  isotropic  material 
(reference  5). 

f(2)  The  potential  of  obtaining  relative  insensitivity  to  crack- 
propagation  due  to  accidental  damage  in  filamentary  arrays  typ¬ 
ical  for  filament- wound  structures  has  been  shown  theor eticallyj 
by  Hedgepeth  (reference  6).  Furthermore,  tests  on  single  fi¬ 
bers  of  many  materials  show  increased  strength  and  resistance 
to  elevated  temperature  (reference  7),  creep,  and  fatigue  as 
^ompared  to  the  properties  of  the  same  materials  in  bulk  form. 

])  Realization  of  the  implied  potential  performance  gains  depends 
upon  proper  utilization  of  fibers  (isotensoid  design)  and  remains 
to  be  genre. rally  demonstrated  in  practice,  f 


I  (3)  The  peculiar,  non-isotropic  character  of  the  filamentary 
textur^-materiajL  may  be  used  to  advantage  in  specific  applica¬ 
tions^*' One  of  particular  interest  is|the  design  of  expandable 
structures  with  variable  geometr vj  for  which  a  multitude  of 
space  applications  are  presently  considered.  A  general  analyt¬ 
ical  treatment  of  the  problem  of  large  deformations  in  elastic 
materials  reinforced  by  inextensible  chords  has  been  presented 
in  reference  8. 


There  exists  a  need  for  a  fundamental  approach^o  the  problem 
of  design  synthesis  of  filamentary  structures.  The ^pr ki entfTeport  /- 
deals  with^a  specific  and  relatively  important  species  of  felame nt y 
spaicture&J  namely  those  produced  by  filament -winding  processes, 
poth  from  the  point  of  view  of  developing  a  basic  understanding,  and  j 
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ffrom  the  point  of  view  of  particular  applications  to  toroidal  pressure 
containers^  ~  f 

rJ?r^ 

This  investigation  was  performed  with  the  financial  assistance 
of  the  National  Aeronautics  and  Space  Administration 
under  Contract  NASr-8. 
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Fiber  helix  parameter 

Bending  stiffness  of  analog  column 

Coefficients  of  the  first  fundamental  form 

Coefficients  of  the  second  fundamental  form 

Elliptical  integral  of  the  second  kind 

Force/unit  length 

Elliptical  integral  of  the  first  kind 

Pressure  load  parameter 

Fraction  of  filament  count 

Load,  (also  used  as  "point  on  surface") 

Non-dimensional  radial  coordinate 

Fiber  force 

V  olume 

Radius  vector 

Non-dimensional  axial  coordinate 
Unit  vector 


Distance  between  poles 

Integer 

Curvature 

Argument  of  elliptical  integral 

Specific  strength 

Number  of  filaments  (fiber  count) 

Mass /unit  length 

Pressure 

Radius 

Arc  length 

Curvilinear  coordinates 
Cartesian  coordinates 
Centrifugal  load  parameter 
Slope  of  meridian 
Helix  angle 

Central  angle  of  fiber  spacing 
Central  angle 
Radius  of  curvature 
Angular  velocity 
Integration  variable 
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THEORY  OF  MONO  TROPIC  MEMBRANES 


Definitions 


Consider  a  thin-walled  shell  made  from  a  material  which  con¬ 
sists  of  an  array  of  filaments  as  shown  in  figure  1.  The  thickness  of 
the  shell  may  be  variable,  but  will  always  be  sufficiently  small  so 
that  the  geometrical  properties  of  .the  shell  can  be  described  by  its 
middle  surface  (membrane). 


Consider  now  a  small  but  finite  domain  of  this  shell.  The  geo¬ 
metrical  arrangement  of  the  filaments  traversing  this  domain  shall 
be  such  that  the  tangents  to  all  filaments  have  a  uniquely  determined 
direction,  described  by  an  angle  between  the  filament  centerlines 
and  a  coordinate  line  inscribed  on  the  membrane  surface.  The  cen¬ 
terlines  of  all  filaments,  thus,  belong  to  a  "one  parametric"  family 
of  curves.  Mathematically,  this  can  be  expressed  by 


P  V,  v) 


(i) 


where  u  and  v  are  surface  coordinates  and  .  is  a  unique  and 

(u,  v) 

continuous  function  of  the  surface  coordinates  u  and  v  .  The  con¬ 
dition  of  uniquely  defined  filament  direction  shall  hold  for  all  regular 
points  of  the  membrane,  but  may  not  hold  for  points  of  singularity 
such  as  poles,  lines  of  cross-over,  etc. 

Assume  that  the  filaments  are  completely  flexible  and  capable  of 
carrying  only  normal  forces  only  in  the  direction  jS  of  their  local 
tangents.  If  these  filaments  are  bonded  together  by  a  matrix  mate¬ 
rial,  consider  this  matrix  to  be  completely  compliant  such  that 
stresses  in  the  matrix  are  negligible.  Thus,  the  only  stresses  that 
can  exist  in  this  shell  are  normal  stresses  in  the  locally  defined  di¬ 
rection  /?  . 

A  structure  of  this  nature  shall  be  called  a  Monotropic  Mem- 
brane,  to  indicate  the  unique  and  predetermined  character  of  stress 
distribution. 


Filament- wound  structures  can  frequently  be  considered  as  com¬ 
binations  of  contiguous  or  intermeshing  monotropic  membranes. 
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This  is  particularly  true  if  the  filaments  are  relatively  rigid  and  the 
matrix  or  binder  material  is  relatively  flexible,  so  that  the  predom¬ 
inant  part  of  the  structural  loading  is  carried  by  the  filamentary 
constituent  of  the  composite  material. 


Geometry  of  Monotropic  Membranes 

In  this  paragraph,  geometrical  relations  referring  to  surfaces 
and  curves  on  surfaces  are  summarized  from  references  9  through 
11  for  convenience  of  the  subsequent  discussion  of  the  properties  of 
monotropic  membranes. 

Consider  the  middle  surface  of  the  shell  described  by  a  vector 
X  extending  from  an  arbitrary  origin  MOn  to  any  point  nPM  on  the 
surface  shown  in  figure  2.  The  surface  is  inscribed  with  a  system 
of  orthogonal  curvilinear  coordinates  u,  v  .  For  the  purpose  of 
this  discussion,  the  lines  u  =  Constant  and  v  =  Constant  will  be 
chosen  as  the  lines  of  principal  curvature  of  the  surface* 

Let  (e  ,"e  ,*e  )  be  a  triad  of  mutually  orthogonal  unit  vectors, 
t  n  g 

moving  along  the  centerline  of  a  filament  located  on  the  surface 
where 


e 


t 


is  tangential  to  the  filament  centerline 


e  is  perpendicular  to  the  surface 

e  is  located  in  the  tangent  plane  to  the  surface 
®  centerline  and  normal  to  the  filament  tangent. 


The  rate  of  change  of  the  unit  tangent  vector  e^  ,  while  pro¬ 
ceeding  along  the  filament  centerline  for  a  differential  length  ds  , 
is  equal  to  the  curvature  vector: 


k 

where  k  is  a  vector  collinear 
and  of  absolute  magnitude  ^  . 


de^ 

ds 


(2) 


with  the  radius  vector  of  curvature  p 
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The  normal  curvature  vector  k  is  the  projection  of  k  into  a 


plane  spanned  by  e  and  e 

t  n 


Its  magnitude,  the  normal  curvature, 
can  be  expressed  by  the  scalar  vector  product 


n 


=  <k'*n> 


(3) 


The  geodesic  curvature  vector  k^ 


is  the  projection  of  k  into 


the  tangent  plane  to  the  surface.  Its  magnitude,  the  geodesic  curva¬ 
ture,  can  be  expressed  by  the  scalar  vector  product 


k  =  (k  •  e  )  (4) 

g  g 


Curves  of  vanishing  geodesic  curvature  are  called  geodesics.  They 
have  the  characteristics  of  straight  lines  in  the  non- Euclidian  ge¬ 
ometry  of  the  surface  described  by  X  .  Thus,  they  are  also  curves 
of  minimum  (or  maximum)  distance  between  two  points  on  the  sur¬ 
face.  A  particularly  important  characteristic  of  geodesics  is  the 
fact  that  geodesics  retain  their  geodesic  character  during  an  inex- 
tensional  (bending)  deformation  of  the  membrane. 

The  pertinent  geometrical  properties  of  the  surface  can  be  ex¬ 
pressed  in  terms  of  two  fundamental  quadratic  forms: 

(1)  The  first  fundamental  form  describes  the  length  ds  of  a 
line  element  on  the  surface.  For  an  orthogonal  surface  coor¬ 
dinate  system  (u,  v)  this  first  fundamental  form  reduces  to 

(ds)2  =  (dX  •  dX)  =  E^u2  +  Gjdv2  (5) 

where 


(2)  The  second  fundamental  form  is  a  measure  of  the  change  of 
the  surface  tangent  plane  or  surface  normal  vector,  i.  e.  ,  a 


measure  for  the  curvature  and  twist  of  the  surface.  For  or¬ 
thogonal  coordinate  systems,  the  second  fundamental  form  re 
duces  to 


-  (dX  •  de 


2  2 
Ejjdu  +  G^dv 


where 


Now  let  k  ,  v  and  k  .  .  be  the  normal  and  geodesic  curva- 
n(v)  g(v) 

ture  of  the  lines  v  =  Constant  ,  and  let  k  .  .  and  k  ,  .  be  the 

n(u)  g(u) 

normal  and  geodesic  curvature  of  the  lines  u  =  Constant  . 

The  normal  curvature  of  the  filament  centerline  forming  an  angle  jS 
with  the  line  v  =  Constant  is  then  given  by  Euler's  Theorem  for  the 
orthogonal  system  of  surface  coordinates  formed  by  the  lines  of 
principal  curvature: 

2  2 

k  lQ.  =  k  cos  /3  +  k  sin  £  (7) 

n(p)  n(v)  n(u) 

The  geodesic  curvature  of  the  filament  centerline  forming  an  an¬ 
gle  P  with  the  line  v  =  Constant  is  given  by  Liouville's  Theorem 
for  any  orthogonal  system  of  surface  coordinates: 


k  /m  =  -T-  +  k  ,  N  cosjS  +  k  sin)3  (l 

g(P)  ds  g(v)  g(u) 

Finally,  the  normal  and  geodesic  curvatures  of  the  coordinate 
lines  are  related  to  the  coefficients  of  the  first  and  second  funda¬ 
mental  forms  of  equations  (5)  and  (6)  as  follows: 


(9) 
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afi 

dv 


2Et/G- 


(10) 


Equilibrium  Condition 


Consider  an  element  of  a  filament  of  length  ds  ,  subject  to  an 
axial  force  vector  T  and  an  external  force  vector /unit  length  F' 
of  arbitrary  direction  as  shown  in  figure  3.  Equilibrium  for  the 
element  demands  that 


T(  s+ds) 


ds 


0 


or,  in  differential  vector  form 
dT 


s) 


ds 


+  F'  =  0 

s 


(ID 


The  vector  T^  can  be  written  as  -  e^  •  ,  where  T^ 

is  the  magnitude  of  the  force  vector  T  at  s  .  Thus,  equation 
(11)  can  be  written  as 


_d 

ds 


et)  +  F' 


0 


and  by  partial  differentiation  of  the  first  term,  considering  equa¬ 
tion  (2)  yields: 


dT/  x 

-r1  + 

ds  t 


T.  .  k  +  F'  .  =  0 

(s)  (s) 


(12) 


Scalar  multiplication  of  equation  (12)  in  turn  with  the  three  unit 

vectors  of  the  triad  (e  .e  ,e  )  yields  three  scalar  equilibrium 

t  n  g 


conditions : 
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where 


dT 

— . —  -f 

F' 

— 

ds 

t 

Tk  + 

F' 

_ 

n 

n 

Tk  + 

F' 

_ 

g 

g 

4 

c-h 

II 

(F'- 

V 

F*  = 

(F'  . 

e  ) 

n 

n 

F'  = 

(F'. 

e) 

g  g 


(13) 

(14) 

(15) 


are  the  three  components  of  the  external  load  vector  in  the  three 
directions  indicated  by  the  triad  (e  ,  e  ,e  )  ,  respectively  (fig¬ 
ure  3).  t  n  g 


Discussion 


The  equations  (13)  through  (15)  relate  the  filament  load  to  the 

geometrical  characteristics  of  the  filament  given  by  e  ,  k  , 

t  n 

and  k  .  The  two  curvatures  are  related  to  the  properties  of  the 


lines  of  principal  curvature  serving  as  surface  coordinates  by 
Euler's  and  Liouville's  Theorems,  equations  (7)  and  (8)  respec¬ 
tively,  and  may  be  related  to  the  coefficients  of  the  fundamental 
forms  as  given  in  equations  (9)  and  (10). 


An  examination  of  the  equilibrium  equations  (13)  through  (15) 
yields  the  following  results: 


(1)  For  a  given  monotropic  shell,  only  one  of  the  three  com¬ 
ponents  of  external  load  vector  F'  can  be  freely  chosen,  while 
still  satisfying  all  conditions  of  equilibrium.  This  implies  that 
the  monotropic  shell  is  twofold  statically  underdeterminate, 
i.  e.  ,  it  will  act  as  a  mechanism  rather  than  as  a  structure  un¬ 
der  all  external  loads  not  satisfying  the  conditions  of  equations 
(13)  through  (15).  This  effect  may  be  beneficial  where  a 


11 


"variable  geometry  structure"  such  as  a  foldable  and  expandable 
shell  is  desired.  On  the  other  hand,  however,  peculiar  insta¬ 
bility  phenomena,  as  discussed  in  the  latter  portion  of  this  re¬ 
port,  may  occur,  requiring  particular  study  and  attention. 

In  general,  it  is  necessary  to  have  three  contiguous  or  in¬ 
termeshing  layers  of  monotropic  membranes  present  in  a  struc¬ 
tural  shell  in  order  to  produce  a  statically  determinate  system. 

(2)  For  the  case  where  F’^  =  0  ,  the  fiber  force  T  is  invari¬ 
ant  throughout  the  structure.  Such  a  structure  is  called  "Iso- 
tensoid"  (from  iso- equal,  tensus- extension)  to  indicate  the 
uniform  character  of  stress  and  strain  in  the  load  carrying  ma¬ 
terial.  An  isotensoid  structure,  furthermore,  constitutes  an 
optimum  design,  since  the  material  will  be  stressed  equally 

(i.  e.  ,  utilized  uniformly)  at  all  points. 

(3)  For  monotropic  shells  that  are  loaded  by  pressure  loads 

normal  to  the  surface,  both  F1  and  F1  vanish.  The  struc- 

t  g 

ture  is  isotensoid  and  the  geodesic  curvature  k  of  the  fila- 

g 

ments  vanishes;  therefore,  the  curves  of  filament  centerlines 
describe  geodesics  of  the  membrane  surface.  These  are  of 
particular  importance  for  the  case  of  filament  winding,  since  a 
fiber  placed  under  some  tension  on  an  arbitrarily  shaped  man¬ 
drel  will  tend  to  follow  a  geodesic  on  the  mandrel  surface. 

The  general  problem  of  optimum  isotensoid  design  synthesis  for 
a  filament-wound  structure  may  now  be  stated  as  follows: 

For  a  given  external  load,  find  the  surface  shape  and 
associated  filamentary  geometry  of  one  or  several  mono¬ 
tropic  membranes  for  which  F'  vanishes. 

This  problem,  applied  to  pressure  loaded  shells  of  revolution, 
will  be  discussed  in  detail  in  the  remaining  portion  of  the  report.  A 
different  situation  arising  from  a  combination  of  pressure  and  cen¬ 
trifugal  loads  is  discussed  in  Appendix  A.  This  case  demonstrates 
the  application  of  the  basic  approach  to  filamentary  design  of  not 
necessarily  isotensoid  character  and  for  more  complex  loading  con¬ 
ditions. 


12 


ISOTENSOID  PRESSURE  SHELLS  OF  REVOLUTION 


Basic  Equations 

Consider  a  surface  of  revolution  (r,  <p,  z)  ,  as  shown  in  figure  4. 
The  coordinates  of  the  vector  X  extending  from  the  origin  of  a  Car¬ 
tesian  system  (x  y  z)  are 


X 


(16) 


The  lines  of  principal  curvature  used  as  surface  coordinate 
lines  are  the  parallel  circles  u  =  z  =  Constant  and  the  meridians 
v  =  <p  =  Constant  . 

The  differential  vector  dX  in  terms  of  the  coordinates  (zcp)  , 
is : 


The  unit  vector  normal  to  the  surface  has  the  coordinates 


e 

n 


(18) 


where  a  is  the  angle  between  the  local  surface  normal  e  and  the 
(x  y)  plane.  n 

The  second  derivatives  of  X  ,  with  respect  to  the  principal 
curvature  coordinates,  are: 
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“  cos<p 


smcp 


-  r  cos  <p 

-  r  sin tp 

0 


From  equations  (5)  and  (17)  the  coefficients  of  the  first  fundamental 
form  become: 


,  /  dr 

1  +  (  dz 


SX 

3X 

dz 

dz 

8X 

ax 

d(p 

b<p 

Similarly,  from  equations  (6),  (18),  and  (19)  the  coefficients  of  the 
second  fundamental  form  become: 


a2x  _ 


82X  _ 

— T  ’  e 


r  cosa 


From  equations  (20)  and  (21)  the  principal  curvatures  of  the  coor¬ 
dinate  lines  become: 
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(£)' 


2  \  3/2 


—  (22) 


cosa 


The  negative  signs  appearing  in  equations  (22)  and  (23)  are  a  con¬ 
sequence  of  the  sign  conventions  used  for  a  and  O,  ,  as  indicated 

M 

in  figure  4.  It  is,  of  course,  possible  to  derive  the  values  of  prin¬ 
cipal  curvature  by  inspection  of  figure  4  directly.  The  analytical 
development  has  been  given  here  to  demonstrate  the  application  of 
an  approach  that  may  be  used  for  analysis  of  generalized  surfaces. 

The  normal  curvature  in  direction  /?  becomes,  according  to 
Euler's  Theorem  (equation  (7)): 


cos  B  ~ 


cosa  .  2- 

-  sin  8 

r 


The  direction  of  the  geodesic  filament  is  obtained  by  integration 

of  Lriouville's  Theorem  (equation  (8))  with  k  Q  =  0  .  For  surfaces 

g(p) 

of  revolution,  this  integration  yields  the  Theorem  of  Clairaut: 


sinfl  =  sin8  —  (2 

or  x 

where  sin^  and  r^  are  arbitrary  constants  of  integration. 

Let  C  =  sin/?  :  filament  helix  parameter  characteristic 

for  the  filamentary  geometry. 


R  =  r/r  :  non-dimensional  radius  of  shell, 
o 


Then,  equation  (25)  can  be  written  as 
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(25a) 

^R2  -  C2  (25b) 

Consider  now  a  structural  shell  of  revolution  consisting 
of  several  layers  of  monotropic  membranes,  as  shown  in  figure  5, 
subject  to  a  rotationally  symmetrical  pressure  load  p^  that  may 

be  variable  along  the  meridian,  but  is  uniform  along  any  parallel 
circle.  Each  monotropic  membrane  is  assumed  to  comprise  n^ 

filaments  uniformly  spaced  along  any  parallel  circle  of  the  shell, 
each  filament  of  each  monotropic  membrane  forming  an  angle 
sinjS  =  C.  with  the  meridian  at  the  extreme  perimeter  of  radius  r  , 
i 

and  each  filament  carrying  an  equal  and  uniform  axial  force  T. 

The  normal  load/unit  length  carried  by  each  filament  is 


sinjS  = 

and 

cos/3  = 


F* 


n. 

l 


2r  77 


COSj8. 


(26) 


where  p.  is  the  normal  pressure  carried  by  each  monotropic 
membrane  "i"  . 


The  equilibrium  condition  (equation  (14)),  considering  equations 
(26)  and  (24),  then  becomes 


2 

,  sin  B.  ~  _ 

1  0  cosa  .  2r7T  „ 

- cosp.  -  -  - 3—  +  p. — —  =  0 

O, ,  1  r  cosfi.  in.  T 

M  1  1 


(27) 


n. 

introducing  the  Mnon-dimensionalM  filament  count  N ^  =  —  ,  where 

n  is  the  total  number  of  filaments  traversing  a  parallel  circle,  and 
substituting  equations  (25a)  and  (25b)  into  equation  (27)  yields: 


pi 


nT 

2 

277r  R 
o 


(28) 


Summing  over  all  monotropic  membranes  "iM 
non-dimensional  meridional  curvature  yields: 


and  solving  for  the 


sna/r2  -  ciZ 


osa  v' 

R  Z 


N.  C. 
1  1 


2  P(z) 


277  p  r 

where  K  = - -  is  the  non-dimensional  pressure  load  param¬ 

eter  which  is  characteristic  for  the  curvature  properties  of  the  shell. 

Equation  (29)  describes  the  meridional  curve  of  the  shell  of  rev¬ 
olution  in  terms  of  the  intrinsic  coordinates  n  and  OL  .  as  a 

M 

function  of  the  characteristic  parameters  C.  and  K  . 

An  integral  form  of  equation  (29)  can  be  derived  by  considering 
the  equilibrium  conditions  pertaining  to  the  whole  shell,  as  shown  in 
figure  5. 


Equilibrium  in  axial  direction  requires: 


|£N.  nT  cos£^  |  -  r^  77pQ  =  |LN.  nT  cosce  coSjS^j  -  R^r^TTp 


or,  solved  for  costt: 


cosa  = 


E  N  , /r2  -  C.2 


sNi\/^7)- 1  \i-r2\ t 


Equilibrium  in  torque  around  the  z-axis  for  a  single  filament 
requires  (since  the  pressure  loading  can  not  contribute  to  a  torque) 


T  sinfl.  r 

l 


T  sin B  r 

o.  o 

i 


which  simply  reduces  to  the  Theorem  of  Clairaut  (equation  (25))  by 
division  through  T. 
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For  the  complete  shell,  torsional  equilibrium  in  absence  of  an 
externally  applied  torque  requires: 

EN.  C.  =  0  (31) 

i  i 

A  further  condition  is  imposed  for  closed  pressure  vessels  of 
singly  connected  volume  without  externally  applied  axial  forces 
("Bottles"). 

For  this  case,  axial  equilibrium  requires  that: 

EN.  yi  -  C.  =  (32) 

l  v  l  2 

Finally,  the  problem  of  defining  the  filamentary  geometry  can 
be  solved  by  integration  of  the  differential  geometrical  relationship 
between  <p  and  a  ,  along  a  geodesic.  By  inspection  of  figure  4 
we  obtain: 


da  tan/?  =  r  d<p 


or,  integrated 

V) 


<P 


(a  =  o) 


(33) 


Methods  of  Solution 


Equation  (30)  represents  an  integral  form  of  equation  (29) 
which  can  be  transformed  into  a  first-order  non-linear  differential 
equation  for  R  suitable  for  numerical  integration  by  substituting 


cosa  = 


2 


where  Z  =  —  is  the  non-dimensionalized  axial  coordinate, 
r 

o 
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The  res vilting  differential  equation  can  be  written  in  finite  differ¬ 
ence  form  for  numerical  integration  by  means  of  a  digital  computer. 


A  program  has  been  prepared  for  an  IBM,  Type  1620  computer, 
capable  of  handling  arbitrary  inputs  for  N.,  C.  and  K  ,  and  com¬ 
puting  both  intrinsic  (p,  a)  and  non-dimensionalized  Cartesian  co¬ 
ordinates  (R,  Z)  of  the  meridional  shapes  by  means  of  a  modified 
Runge-Kutta  method  (reference  12).  Additional  output  data  available 
from  the  program  are  the  filament  geometry  <p  as  well  as  sur- 


faces  S  and  enclosed  volumes  Y  of  the  bodies  of  revolution  gen¬ 
erated  by  rotating  the  meridional  curves  (R,  Z)  around  the  z-axis. 


Analytical  solutions  of  equation  (30)  are  possible.  These  lead  to 
expressions  which  may  be  reduced  to  elliptical  integrals.  In  the 
general  case,  however,  the  resulting  expressions  are  quite  involved 
and  time-consuming  for  evaluation.  As  a  consequence,  direct  digital 
integration  appears  to  be  a  more  economical  approach  for  the  pur¬ 
pose  of  obtaining  design  data.  An  exception  is  the  case  for  C  =  0  , 
discussed  in  Appendix  B.  There,  the  solution  can  be  reduced  to 
simple  combinations  of  elliptical  integrals  of  the  first  and  second 
kind,  which  are  tabulated,  for  instance,  in  reference  13. 


For  the  purpose  of  a  graphical  solution,  equations  (29)  and  (30) 

are  suitable  in  their  intrinsic  form.  Tables  of  p  and  a  in  func- 

M 

tion  of  r/r  can  be  computed  for  any  given  set  of  K,  C  ,  N  from 
°  i  i 

equations  (30)  and  (29).  With  slope  0 !  and  radius  of  curvature  p 

M 

known,  very  accurate  meridional  shapes  can  be  constructed  by 
means  of  rule,  compass,  and  protractor. 

In  addition  to  the  numerical  and  graphical  methods  mentioned, 
the  particular  case  of  C  =  0  (i.  e.  ,  meridional  fiber  direction  only) 
allows  solution  of  the  isotensoid  meridional  shape  equations  by 
means  of  a  column  analog,  as  shown  in  figures  6  and  7.  This  analog 
is  described  in  more  detail  in  Appendix  B. 


All  three  methods,  digital  integration,  graphical  solution,  and 
column  analogy,  have  been  used  to  generate  the  data  shown  in  fig¬ 
ures  10  through  20  of  this  Technical  Note. 


Discussion  of  Results 


For  the  purpose  of  an  initial  evaluation  of  possible  isotensoid 
pressure  vessels,  the  number  of  free  parameters  may  be  reduced 
as  follows : 


(1)  The  pressure  p  is  assumed  to  be  constant  throughout 
the  length  of  the  shell,  i.  e.  ,  ^  =  PQ  -  Constant. 


(2)  The  shells  are  considered  to  consist  of  two  monotropic 
membranes  satisfying  the  torque  equilibrium  equation  (31): 


where  =  N  =  0.  5 

and  Ci  =  -  =  C 


With  these  assumptions,  the  intrinsic  equations  (29)  and  (30) 
reduce  to 


(29b) 


(3  0b) 


where  C  and  K  are  the  defining  parameters  of  the  isotensoid 
pressure  shell. 


Inspection  of  equations  (29b)  and  (30b)  shows  that  no  real  so¬ 
lutions  for  either  p  or  Ot  exist  if  R  <  C  .  Thus,  the  minimum 
distance  from  the  axis  of  rotation  which  a  geodesic  filament  can 
have  is  given  by  its  helix  angle  defined  by  C  .  An  important  de¬ 
sign  implication  that  may  be  derived  from  this  is  that  there  will 
always  be  a  polar  opening  for  any  geodesically  helix  wound  (C  /  0) 
pressure  vessel.  The  size  of  the  polar  opening  R  can  be  com¬ 

puted  from  equation  (30b)  by  letting 
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Excluding  the  trivial  solution  R  =  1  ,  a  relation  between 

R  .  is  obtained: 
mm 


K 

2 


K, 


C,  and 


(34) 


This  equation  is  useful  if  R^  given  as  a  design  parameter, 

and  a  compatible  pair  C  and  K  need  to  be  determined  for  the  de¬ 
sign  of  the  pressure  vessel. 


R 


(Z) 


The  solutions  that  may  be  obtained  for  the  meridional  shapes 
by  solution  of  the  intrinsic  equations  (29b)  and  (30b)  can  be 


classified  according  to  their  general  character  for  a  range  of 
parameter  K  . 


For  the  purpose  of  this  discussion,  only  the  range  of  K  ^  0 
will  be  considered.  Negative  K- values  yield  similar  shapes  (cor¬ 
responding  to  either  negative  pressure  or  compressive  fiber  force) 
and  yield  no  new  general  insight  into  the  problem  even  though  they 
may  be  important  for  particular  design  objectives.  The  various 
shapes  for  positive  K- values  are  discussed  below  and  are  shown  in 
schematic  form  in  figures  8  and  9: 


(1)  Hyperboloid  (figure  8(a)) 

Here  the  extreme  perimeter  r  is  a  minimum.  For  all 

positive  z  ,  cos#  is  always  negative  and  is  always 

positive.  This  requires  K  =  0  and  results  in  the  fila¬ 
ments  being  straight  lines  rotating  around  the  z-axis  , 
forming  a  hyperboloid  of  revolution.  This  is  a  surface  of 
zero  curvature,  corresponding  to  zero  internal  pressure  p 


(2)  Cylinder  (figure  8(b)) 


For  a  cylinder  R  =  1  ,  cos#  =  0  and  p^.  =  03  .  These 

conditions  are  satisfied  if 


v>  -  c2 


K 
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(3)  Corrugated  Tube  (figure  8(c)) 

A  corrugated  tube  will  exhibit  a  point  of  inflection  where 
p^.  =  ±  00  at  a  value  of  R  >  C  .  a.  will  always  be  smaller 

7 T 

than  ±  —  .  This  requires  that 

4 


0  <  K  < 


1  -  C 


If  K  < 


,  then  the  corrugations  are  located  out- 


1  -  C 


side  the  cylinder  R  =  1  and  r  is  a  minimum;  if 

c2  ° 

K  > —  ■  ,  then  the  corrugations  are  inside  the 

Vl  -  C2 

cylinder  R  =  1  ,  and  r  is  a  maximum,  as  indicated  in 

figure  8(c). 

For  those  corrugated  tubes  ("bottles”)  which  satisfy 
the  axial  equilibrium  equation  (32),  the  following  relations 
hold: 


The  minimum  radius  R  .  is  related  to  the  fiber 

mm 


helix  parameter  C  by 


1  -  C 


R  .  2  -  C2 

mm 

R  .  3 
mm 


The  slope  O'  of  the  meridian  is  given  by 


cosQf  =  R 


3V^ 

v  R  -  C 


The  radius  at  the  inflection  point  of  the  meridian  is 


R.  ,  =  yfT/2  C 

(P  =  °°) 


22 


The  slope  at  the  inflection  point  is  given  by 

costW)  =  !  s*  c2V‘  - c2 

(4)  Cusp  (figure  8(d)) 

For  the  case  of  a  cusp,  the  angle  a  becomes  indefinite 
an<^  becomes  0  at  R  =  C  .  This  condition  is  satis¬ 

fied  if 


V1  -  c2 

The  cusp  shape  constitutes  the  transition  between  the  cor¬ 
rugated  tube  and  the  progressive  loop. 

(5)  Loop  (figure  9) 

Looped  curves,  which  are  of  particular  interest  for  the 
generation  of  toroidal  pressure  vessels,  are  generated  if 
PM  is  negative  for  all  values  of  R  .  This  requires 


V1  -  C2 

Three  distinct  types  of  loops  are  generated  as  K  increases 

(a)  A  progressive,  periodic  loop  (figure  9(a)). 

(b)  A  closed  loop  forming  the  transition  between  pro¬ 
gressive  and  regressive,  periodic  loops  (figure  9(b)). 

(c)  A  regressive,  periodic  loop  (figure  9(c)). 

The  radius  R  at  which  the  meridional  curve  is  horizontal 
ri 

(i.  e. ,  the  location  of  maximum  thickness  of  the  toroid  gen¬ 
erated  by  rotation  of  the  loop)  is 

R(a=TT/ 2)  =  ^  ’  K  V1  "  C 
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A  summary  of  the  domains  in  the  positive  K- C  plane  corre¬ 
sponding  to  the  various  types  of  solutions  is  shown  in  figure  10. 

Also  shown  in  figure  10  are  the  curves  R  .  =  Constant  ,  and  the 

mm 

C-K  values  resulting  in  MbottlesM  according  to  equation  (32). 

A  number  of  cases  for  various  parameter  combinations  have 
been  worked  out  by  digital  integration  of  equations  (30)  and  (33)  com¬ 
plemented  by  graphical  and  column  analog  studies.  Primary  atten¬ 
tion  was  given  to  the  looped  type  of  solutions  which  are  of  interest  in 
the  design  of  toroidal  pressure  vessels. 

Figures  11  through  14  show  the  meridional  shapes  for  1/2  pe¬ 
riod  (0  ~  a  =  7T)  of  double  helix  wound  shells  according  to  equations 

2 

(29b)  and  (30b)  for  K  > -  -  and  for  C-values  of  0,  .25,  .  5, 

V1  -  °2 

and  .75,  respectively. 

Figures  15  through  17  show  solutions  for  the  case  of  looped 
meridians  for  shells  consisting  of  three  layers;  =  0  ,  =  .  5  , 

and  =  -.  5  ,  with  variable  ratios  of  N^/N^  and  with  =  N3  . 

Figure  18  shows  a  case  with  6  layers,  representing  the  case  of 
a  corrugated  tube,  together  with  the  filamentary  geometry  of  the 
various  layers. 

In  the  cases  of  multiple  C-values,  shown  in  figures  15  through 
18,  return  points  occur  at  R  =  C^  ,  where  the  corresponding  mono¬ 
tropic  membrane  folds  back  into  itself,  reversing  the  helix  angle 
and  thus  forming  C.  ,  ,  =  -  C.  . 

Figure  19  shows  an  example  of  an  unsymmetrical  loop  resulting 
from  linearly  variable  pressure  simulating  the  situation  of  a  pres¬ 
sure  vessel  subject  to  hydrostatic  pressure  gradients  such  as  may 
be  encountered  in  a  liquid  container  during  the  boost  phase  of  a 
rocket. 

Figure  20  shows  a  realization  of  the  particular  case  for  a  closed 
loop  indicating  the  general  shape  and  the  disposition  of  the  winding 
pattern. 
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MORPHOLOGY  OF  CLOSED  PRESSURE  VESSELS 


The  range  of  possible  isotensoid  pressure  shells  of  revolution 
with  meridional  shapes,  as  shown  in  figures  11  through  20,  is  lim¬ 
ited  by  the  restriction  of  the  parameter  ranges  and  combinations  that 
have  been  considered.  Nevertheless,  five  basic  types  of  possible 
closed  pressure  vessel  designs  are  identifiable  from  a  review  of  the 
presented  material.  A  brief  discussion  of  their  general  character¬ 
istics  follows : 


Containers  with  Simply  Connected  Volume 

(1)  Type  1  (figure  21):  Closed  shells  containing  a  monotropic 
membrane  with  meridional  fiber  direction  (C  =  0)  cross¬ 
ing  the  z-axis  . 

In  the  special  case  where  the  "bottle11  relationship  of  K 
and  CL  given  in  equation  (32)  is  satisfied,  the  meridian 

crosses  the  z-axis  at  right  angles,  and  simple  closed 
"bottles"  result.  Where  the  parameters  K  and  CL  do  not 

satisfy  equation  (32),  additional  elements  in  the  form  of  an 
axial  compression  column  or  tension  rod  located  in  the 
z-axis  are  necessary  for  equilibrium.  A  discussion  of 
such  shapes  has  been  given  in  reference  1. 

Shells  containing  monotropic  membranes  with  C  =  0  are, 

from  a  practical  point  of  view,  not  realizable  in  pure  form 

since  all  fibers  of  the  C  =  0  layers  cross  the  z-axis  at 

the  same  location,  resulting  in  a  "pole"  singularity,  and 

generating  an  impractical  local  buildup  of  material  in  the 

winding  pattern  at  the  pole.  This  form  of  pressure  vessel 

has  been  approached  by  employing  very  small  helix  angles 

(or  small  C-values).  Such  structures  have  been  made 
o 

in  the  form  of  the  "Bermuda  bottles"  ( C  ~  0,  K  ~  2)  ,  and 

as  endclosures  of  " end- over- end"  wound  cylinders  pro¬ 
duced  by  several  manufacturers  of  filament- wound  pressure 
vessels. 
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Other  possibilities  in  reducing  the  polar  problem  rest  in  the 
use  of  wide  bands  for  the  C  =  0  layers. 

(2)  Type  II  (figure  22):  Closed  bottles  of  singly  connected  vol¬ 
ume  incorporating  a  secondary  closure  structure  for  the 
polar  openings. 

This  type  of  pressure  vessel  needs  to  satisfy  equation  (32) 
for  axial  balance  unless  axial  tension  or  compression  mem¬ 
bers  supporting  unbalanced  endclosure  forces  are  added  to 
the  design.  As  seen  from  figure  10,  all  "bottles11  fall  in  the 
realm  of  monotropic  membranes  in  the  form  of  corrugated 
tubes.  These  corrugated  tubes  can  be  modified  to  include 
a  cylindrical  portion  of  R  =  1  between  two  converging  por¬ 
tions  by  application  of  a  circumferential  winding  pattern 
(C  =  1)  .  Such  "bottles"  represent  the  ideal  design  for  the 

familiar  type  of  filament- wound  cylindrical  pressure  ves¬ 
sels  and  rocket  motor  cases  with  "ovaloid  endclosures" 
(reference  2). 

A  section  of  a  "corrugated  tube"  will  have  two  polar  open¬ 
ings  with  a  minimum  radius  of  R  .  >  C.  .  If  a  polar 

&  mm  i 

mm 

opening  R  >  C.  is  required  (for  instance,  for  nozzle 

min 

attachments,  etc.  ),  then  a  circumferential  retaining  hoop 
(C  =  1)  must  be  inserted  into  the  winding  pattern  to  obtain 
balance  of  an  isotensoid  structure  requiring  a  special  wind¬ 
ing  arrangement.  In  either  case,  a  secondary  closure 
structure,  normally  consisting  of  a  flanged  ring  or  disk  in¬ 
sert,  must  be  provided. 

A  discussion  of  the  detail  design  problems  for  the  flange  at¬ 
tachments  is  beyond  the  scope  of  this  report.  Suffice  it  to 
note  that  the  proper  flange  shapes  may  be  derived  from  the 
basic  equations  for  isotensoid  shells  by  introducing  the  con¬ 
tact  pressure  between  flange  and  filamentary  shell  into  the 
equations . 

As  in  the  case  for  Type  I,  excessive  buildup  of  filaments  at 
the  polar  opening  (i.  e.  ,  at  the  "turn  around"  point  of  the 
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winding  pattern)  may  result.  This  problem  can  be  reduced 
by  use  of  a  tension  hoop,  as  discussed  above,  which  pre¬ 
vents  the  filaments  from  piling  up  at  the  turn  around.  Other 
solutions  that  have  found  practical  application  consist  of 
using  filaments  of  finite  width  (ribbon- winding)  or  by  small 
variations  of  the  filament  helix  angle  around  a  mean  value 
(wobble-winding),  which  tends  to  spread  the  buildup  over  a 
finite  domain  of  the  interface  area. 


Containers  with  Doubly  Connected  Volume 

TyPe  (figure  23):  Closed  toroids  (doubly  connected  vol¬ 
umes)  generated  from  progressive  loops. 

Progressive  loops,  as  shown  in  figure  9(a),  may  be  used  to 
generate  close  toroidal  pressure  vessels  in  two  ways: 

(a)  A  compression-resistant  hub  structure  may  be  used  to 

provide  internal  closure  of  a  loop  extending  between 

two  minimum  perimeters  R  .  of  the  meridional 

mm 

curve.  The  filaments,  in  this  case,  may  either  be  con¬ 
tinuous  around  the  torus  meridian,  or  they  may  return 
around  retaining  hoops  located  at  the  hub  rims  in  the 
fashion  of  an  automobile  tire,  similar  to  the  one  shown 
for  Type  II  in  figure  22. 

(b)  A  compression-resistant  hoop  may  be  introduced  in  the 
outward  cusp  of  a  single  loop,  providing  the  necessary 
circumferential  balance  forces. 

In  either  case,  compression- resistant  structural  compo¬ 
nents  are  required  to  close  the  pressure  vessel.  As  a  re¬ 
sult,  the  elastic  deformation  of  the  shell  due  to  pressuri¬ 
zation  cannot  yield  a  geometrically  similar  shape  which  in¬ 
sures  uniform  strain  in  each  filament.  Mismatched  defor¬ 
mations  in  the  various  components  of  the  shell  will  normally 
take  place,  and  deviations  from  ideal  isotensoid  conditions 
will  be  encountered  in  most  cases. 

(2)  Type  IV  (figure  24):  Closed  toroids  from  closed  loops. 

Closed  loops,  as  shown  in  figure  9(b),  form  an  ideal  basis 
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for  the  design  of  isotensoid  toroidal  pressure  vessels,  un¬ 
less  large  access  holes  or  valve  openings  are  required  by 
the  particular  design  application. 

(3)  Type  V  (figure  25):  Closed  toroids  generated  from  regres¬ 
sive  loops. 

Regressive  loops,  as  shown  in  figure  9(c),  may  be  used  to 
form  closed  toroidal  pressure  vessels  in  a  similar  fashion, 
as  discussed  for  Type  III: 

(a)  Tension-resistant  outer  equatorial  bands  of  finite  width 
(C  =  1)  may  be  used  to  provide  external  closure  of  a 
loop  extending  between  two  maximum  perimeters  r 

of  the  meridional  curve. 

(b)  A  tension-resistant  hoop  may  be  introduced  in  the  in¬ 
ward  cusp  of  a  single  loop,  providing  the  necessary 
circumferential  balance  forces. 

(c)  Of  practical  importance  is  the  case  where  both  outer 
bands  and  an  inner  hoop  are  used  for  a  complete  torus 
structure.  The  free  choice  of  width  and/or  spacing  in 
the  outer  bands  allows  greater  freedom  in  obtaining 
particular  design  characteristics.  Further,  the  use  of 
several  bands  allows  construction  of  stable  toroidal 
pressure  vessel  configurations  with  only  meridional 
(C  =  0)  and  perimetral  (C  =  1)  windings. 

As  opposed  to  the  Type  III,  Type  V  employs  tension- 
resistant  material  only,  thus,  the  conditions  of  uniform 
extension  throughout  the  structure  may  be  satisfied  by 
proper  dimensioning  of  the  bands  and  hoops. 

A  review  of  the  various  isotensoid  pressure  vessels  forming 
shells  of  revolution  indicates  that  in  practice  only  the  toroidal  Types 
IV  and  V  can  be  made  as  completely  closed  containers  and  requiring 
only  tension- resistant  filamentary  material. 

These  types  promise  to  be  of  interest  in  two  aspects: 

(1)  Where  foldable,  pneumatically  stabilized  structures  are  re¬ 
quired,  filamentary  materials  made  from  strong  fibers 
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imbedded  in  an  elastic  matrix  may  be  used.  The  absence  of 
any  rigid,  compression-resistant  members  favors  the  Types 
IV  and  V  for  such  application. 

(2)  Secondly,  the  monolithic  structure  of  Types  IV  and  V,  and 
absence  of  additional  weights  due  to  endclosures,  such  as 
required  for  Type  II,  promises  a  high  structural  efficiency 
for  the  toroidal  shapes. 

Toroidal  isotensoid  pressure  vessels  which  were  employed  for 
experimental  verification  of  the  theoretically  derived  data  are  dis¬ 
cussed  below. 


EXPERIMENTAL  STUDIES 


Objectives 

The  experimental  work  reported  here  has  been  accomplished 
with  the  following  objectives  in  mind: 

1.  To  establish  experimentally  the  unrestrained  meridional 
shapes  of  filament- wound  pressure  vessels  and  to  verify 
that  these  shapes  can  be  predicted  analytically. 

2.  To  demonstrate  the  theoretically  predicted  load- deformation 
characteristics  of  toroidal  isotensoid  pressure  vessels. 

3.  To  explore  failure  mechanisms  for  this  type  of  structure. 


Experimental  Method 


General  Approach 

The  objectives  listed  above  have  been  satisfied  by  fabricating  a 
succession  of  filament- wound  toroids,  and  by  measuring  these  under 
pressure.  To  allow  the  shells  to  assume  their  natural,  unrestrained 
shapes,  they  have  been  wound  on  flexible,  inflatable  mandrels  with  a 
fine,  smooth  fiber.  Where  a  bonding  matrix  was  used,  it  has  been 
applied  after  winding  in  the  form  of  liquid  elastomer. 
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The  deformation- vs-load  characteristics  have  been  observed  by- 
subjecting  the  toroids  to  a  sequence  of  internal  pressure  levels. 

Three  distinct  types  of  deformation  have  been  observed: 

(1)  Linear  strains  in  the  fiber  systems. 

(2)  Cross-sectional  shape  changes. 

(3)  Gross  changes  in  the  toroid  shape. 

Material  Selection  and  Design  of  Test  Articles 

For  the  structural  filament,  a  polyester  fiber  in  the  form  of  an 
1100  denier,  untwisted,  250-end  roving  was  chosen.  This  fiber  has 
a  measured  breaking  strength  of  about  15  pounds,  an  ultimate  elonga¬ 
tion  of  about  10  percent,  and  a  cross-sectional  area  equivalent  to 
about  0.  013  inch  solid  round. 

The  matrix  material  used  is  a  polyurethane  elastomer. 

The  reasons  for  the  selection  of  these  materials  are  as  follows: 

(1)  The  high  elongation  under  load  allows  accurate  strain  meas¬ 
urements  to  be  made  by  simple,  direct  means. 

(2)  The  excellent  bending  compliance  of  the  composite  material 
allows  the  pressurized  shell  to  assume  its  ’’natural11  shape, 
unaffected  by  bending  stress  systems  resulting  from  incor¬ 
rect  initial  shape. 

(3)  The  combination  of  high  strength-to- weight  ratio  and  excep¬ 
tional  bending  compliance  of  the  composite  material  used  is 
attractive  for  use  in  efficient  expandable  structures. 

The  major  portion  of  the  test  effort  has  been  devoted  to  studies 
of  the  meridionally- wound  toroid  of  Type  V.  The  size  of  these  to¬ 
roids  was  taken  to  be  that  of  standard  automotive  inner  tubes  since 
this  provides  a  ready  source  of  inflatable  winding  mandrels. 

A  typical  cross-section  for  a  meridionally- wound  toroid  is  shown 
in  figure  26. 
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The  structure  consists  of  three  separate  fiber  systems: 

(1)  A  meridional  winding. 

(2)  An  inner  equatorial  band. 

(3)  An  outer  equatorial  band. 

Other  layers  of  material  shown  in  figure  26  are  non- structural 
and  result  from  the  fabrication  process. 

Design  data  are  given  in  Table  I  for  two  typical  meridionally- 
wound  toroids. 

Fabrication  of  Test  Articles 

The  fabrication  of  the  various  toroids  has  been  accomplished 
with  the  aid  of  the  toroid  winder  shown  in  figure  27.  A  typical  pro¬ 
cedure  is  as  follows: 

(1)  An  automotive  inner  tube  has  been  used  as  a  winding  man¬ 
drel.  The  tube  is  prepared  by  replacing  the  standard  valve 
with  a  flush- mounted  basketball  valve  installed  on  the  outer 
equator.  The  tube  is  then  brought  to  the  correct  meridional 
perimeter  by  wrapping  with  polyester  film  tape  of  controlled 
length  per  turn. 

(2)  The  inner  equatorial  band,  prepared  on  a  separate  mandrel, 
is  then  taped  in  place  on  the  inner  equator  of  the  tube. 

(3)  The  tube  is  covered  with  a  layer  of  light-weight  nylon  cloth 
to  insure  against  premature  failure  due  to  meridional  splits. 

(4)  The  meridional  winding  of  structural  fiber  is  then  applied 
by  the  toroid  winder.  The  fiber  is  either  pre-impregnated 
or  the  winding  is  soaked  with  dilute  elastomer  and  cured  to 
form  the  rubber-fiber  composite  material. 

(5)  The  outer  band,  prepared  separately,  is  then  installed.  The 
valve  hole  in  the  outer  equator  is  kept  clear  by  making  slits 
in  the  windings  parallel  to  the  fibers. 
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Test  Instrumentation 

The  instrumentation  used  in  the  strain- vs-pressure  tests  has 
been  designed  to  measure  the  toroid  deformation  in  three  different 
modes : 

(1)  Linear  elongation  of  typical  fibers  in  each  fiber  system. 

(2)  Volumetric  expansion. 

(3)  Change  of  meridional  shape. 

The  particulars  of  the  various  measurements  are  given  in  Table  II. 

Figure  28  shows  a  diagram  of  the  toroid  cross-section,  indi¬ 
cating  the  various  measurements  which  were  made.  Figure  29  shows 
the  test  set-up  for  toroid  measurements  with  the  extensometer s  used 
to  measure  changes  in  meridional  and  equatorial  lengths  and  the  polar- 
coordinate  ring  used  to  determine  the  meridional  shape. 

The  toroid  is  pressurized  by  a  motor-driven  gear  pump,  which 
transfers  water  from  a  weighable  reservoir  to  the  toroid.  The  volume 
change  is  determined  from  the  reservoir  weight.  The  initial  volume  is 
determined  by  weighing  the  toroid  before  and  after  filling  with  water. 

Test  Procedure 

The  polyester  fiber  used  in  this  test  program  shows  peculiar 
elastic  effects  involving  hysteresis  and  strain  rate  transients.  These 
were  determined  by  a  series  of  load- elongation  tests  made  with  single 
fibers,  bands,  and  cylinders.  Typical  results  from  these  tests  are 
shown  in  figures  30  and  31.  The  effects  of  this  behavior  have  made  it 
necessary  to  standardize  on  a  fixed  schedule  for  tension  vs  time  in  all 
tests  used  for  comparison  of  toroid  tests  to  nominal  materials  prop¬ 
erties  in  order  to  obtain  a  reproducible  tension- vs- strain  diagram. 

The  time  intervals  in  the  tension- vs-time  schedule  for  the  toroid 
tests  were  determined  by  the  requirements  of  the  test  itself.  It  was 
found  that  five  minutes  were  required  to  complete  the  readings  at  each 
pressure  level,  and  that  a  settling  time  of  ten  minutes  was  required 
for  the  strain  to  settle  to  an  essentially  constant  value  over  the  read¬ 
ing  period. 
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In  order  to  determine  the  pressure  which  corresponds  to  a  given 
tension  level,  considering  the  large  deformations  that  occur,  the 
procedure  described  in  Appendix  C  has  been  used,  along  with  a  pre¬ 
diction  of  the  strain  at  each  level  of  tension,  based  on  preliminary 
load- strain  tests  of  the  fiber  composite  material. 


Results 


Instability  Characteristics 

Instability  phenomena  have  been  observed  in  toroidal  pressure 
vessels  and  reported,  for  instance,  in  reference  14.  Such  instabili¬ 
ties  will  be  expected  particularly  in  those  cases  which  are  of  struc¬ 
turally  underdeterminate  nature,  i.  e.  ,  in  shells  containing  two  or 
less  fiber  directions  in  certain  domains  of  the  structure.  The  in¬ 
stability  mechanisms  are  similar  to  those  observed  in  magnetic  bot¬ 
tles  designed  to  confine  plasma  at  high  pressure,  where  the  magnetic 
lines  of  force  take  the  place  of  the  filamentary  texture  (See,  for  in¬ 
stance,  reference  15).  Other  mechanisms  of  this  type  are  those 
found  in  soap  bubbles  (reference  16)  or  may  be  of  the  nature  of  the 
"hydrostatic  instabilities"  observed  in  straight  cylindrical  tubes  made 
from  isotropic  material  (reference  17). 

Early  tests  of  meridionally- wound  toroids  showed  that  similar 
instabilities  are  apt  to  exist  in  these  structures.  Large- deformation 
failure  modes  arise  from  the  tendency  of  the  toroid  to  increase  vol¬ 
ume  at  the  expense  of  locally  high  shear  and  normal  strains  in  the 
matrix  material  which  bonds  the  meridional  fibers  together. 

Two  different  modes  of  this  type  of  buckling  have  been  observed: 

(1)  "Rolling"  Modes 

If  the  toroid  consists  of  a  meridional  winding  which  is  re¬ 
strained  circumferentially  by  an  outer  band  only,  (Type  V(a)), 
the  toroid  will  exhibit  a  tendency  for  the  meridional  loops  to 
roll  over  and  out  of  the  equatorial  band.  This  type  of  defor¬ 
mation  is  shown  in  figure  32. 

(2)  "Baseball  Seam"  Modes 

If  the  toroid  consists  of  a  meridional  winding  restrained  by 
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an  inner  band  only,  (Type  V(b)),  the  toroid  will  buckle  by 
curling  into  a  ’’baseball- seam"  curve.  Several  types  of 
this  mode  are  shown  in  figures  33  through  36,  some  exhib¬ 
iting  striking  geometrical  symmetry  and  beauty. 

Strain  Distribution 

The  results  of  the  strain  measurements  made  during  the  test  of 
Toroid  #4  are  shown  in  figure  37.  The  solid  curve  in  the  figure  is  a 
plot  made  from  tests  of  the  material.  The  plotted  points  represent 
measured  values  of  strain  vs  pressure  (reduced  to  fiber  tension  by 
the  method  outlined  in  Appendix  C) . 

Note  that  the  volumetric  strain  (as  defined  in  Appendix  C)  is  con¬ 
sistently  greater  than  the  strain  in  the  meridian  or  that  in  the  outer 
equator  measured  along  its  average  circumference.  This  observa¬ 
tion  points  to  an  inherent  inadequacy  in  this  particular  design;  with¬ 
out  anomalous  volume  growth,  the  rolling  could  be  attributable,  for 
example,  to  a  combination  of  small  manufacturing  errors  and  high 
sensitivity  of  the  design  to  such  errors. 

The  strain  measurements  obtained  during  test  of  Toroid  #5  are 
shown  in  figure  38  along  with  a  curve  of  the  nominal  properties  of 
the  fiber  obtained  from  tests  under  equal  load  sequence. 

The  slight  tendency  to  roll  is  exhibited  by  the  three  strain  meas¬ 
urements  made  on  the  outer  equatorial  band.  This  roll  distortion 
caused  the  inner  band  to  be  shifted  about  0.  2  inch  from  the  predicted 
plane  of  symmetry  prior  to  failure.  However,  the  volumetric  strain 
was  not  significantly  different  from  the  other  strain  measurements. 

Cross-Sectional  Shapes 

A  comparison  of  the  cross-sections  of  Toroid  #5  at  two  pressure 
levels  with  the  theoretical  curves  is  shown  in  figure  39.  The  heavy 
curve  shown  is  the  theoretical  isotensoid  meridional  torus  section, 
which  is  a  composite  of  two  isotensoid  curves  connected  at  the  out¬ 
side  by  the  outer  band  and  at  the  inside  by  a  rounded  corner  repre¬ 
senting  the  inner  band  contour.  The  light  curve  is  a  geometrically 
similar  enlargement  of  the  heavy  curve,  using  the  tor oidal- coordinate 
origin  as  the  invariant  point.  The  enlargement  factor  is  1.  07  (i.  e.  , 

7  percent  strain)  corresponding  to  a  fiber  tension  of  9  pounds. 
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The  points  plotted  along  these  meridional  curves  represent 
measurements  of  the  cross-section  made  with  the  polar  coordinate 
ring  shown  in  figure  29.  The  points  marked  with  circles  are  from  a 
measurement  at  15  psi  internal  pressure  or  about  0.  5  pounds  per  fi¬ 
ber,  while  the  solid  points  represent  the  cross-section  at  240  psi 
internal  pressure  or  9  pounds  per  fiber. 

In  each  case,  the  measured  points  have  been  located  radially 
with  respect  to  a  point  obtained  from  a  measurement  of  the  outer 
equator  length.  Hence,  the  corresponding  sets  of  experimental  and 
theoretical  cross-sectional  data  have  been  placed  on  the  diagram  in¬ 
dependently  from  each  other;  these  data  are  related  through  the 
curve  of  nominal  fiber  properties  (i.  e.  ,  7  percent  strain  at  9  pounds 
tension)  but  otherwise  bear  no  direct  relation. 

Examination  of  the  data  obtained  at  low  pressure  reveals  that 
the  meridian  is  initially  more  nearly  circular  than  theoretically  pre¬ 
dicted.  This  is  attributable  to  the  fabrication  process :  the  matrix 
material  was  applied  and  allowed  to  cure  at  low  internal  pressure  of 
the  toroid.  The  initial  shape,  therefore,  is  influenced  by  the  inner- 
tube  which  tends  to  approach  a  circular  section. 

At  the  higher  pressure  level,  the  load  carried  by  the  fiber  is 
large  compared  with  the  load  carried  transversely  by  the  matrix  and 
rubber  tube,  and  the  .meridional  shape  is  determined  primarily  by 
the  fiber.  In  this  case  the  comparison  between  the  theoretically  pre¬ 
dicted  and  the  measured  meridional  shapes  shows  excellentragree- 
ment.  In  the  full-scale  layout,  the  experimental  points  all  fall  within 
.  03  inch  of  the  theoretical  curve. 

Additional  cross-section  measurements  have  been  made  with  a 
toroid  of  Type  V(a)  with  inner  band  only  (Toroid  #3),  shown  com¬ 
pleted  in  figure  40,  and  with  a  helix- wound  toroid  of  Type  IV  (To¬ 
roid  #6)  shown  in  the  process  of  fabrication  in  figure  27. 

Fair  agreement  between  measured  and  theoretical  (MisotensoidM) 
meridional  cross-sections  for  both  toroids  are  shown  in  figures  41 
and  42. 

Astro  Research  Corporation 

Santa  Barbara,  California,  July  20,  1961 
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APPENDIX  A 

MONOTROPIC  SHELLS  SUBJECT  TO 
CENTRIFUGAL  LOADING  AND  PRESSURE 


Consider  a  monotropic  membrane  in  the  form  of  a  shell  of  revo¬ 
lution  similar  to  the  one  shown  in  figure  5. 


In  addition  to  the  normal  load / unit  length  of  filament  due  to  a 

uniform  internal  pressure,  also  consider  a  radial  load  vector  F1 

to 

caused  by  a  rotation  of  the  shell  around  the  z-axis  with  an  angular 
velocity  CO  .  /  \ 

(cost p) 


F'  =  m'  r 
CO 


CO  /sin<p\ 

\  0  ) 


(1-A) 


where  m1  is  the  mass  associated  with  a  unit  length  of  the  filament. 


The  tangential  component  of  this  load  vector  is  obtained  by  sea 


lar  multiplication  with  the  unit  tangent  vector  e^  : 


*'  =  (f'  •  e  ) 
t  V  co  t/ 


=  m'  rCO  cos0  sinO? 


(2- A) 


The  first  equilibrium  equation  (13)  becomes 


+  m'  rCO  cos/3  sintt  =  0 


and,  since  dr  =  -  ds  cos/3  sintt  : 


dT  =  m'  rco  dr 


which  integrates  to 


T  -  T 


2  1/2  2  \ 

m  “  2  (ro  -  r  ) 


(>  - 


(3-A) 
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m'  00  r 

where  T  is  the  fiber  force  at  r  and  =  — -= -  is  the  non- 

o  o  T 

o 

dimensionalized  centrifugal  load  parameter.  Thus,  the  fiber  force 
is  non-uniform  and  an  nisotensoidM  design  would  require  a  tapering 
cross-section  of  filaments  designed  to  keep  the  filament  stress  and 
strain  uniform. 

The  equilibrium  conditions  for  the  shell  of  revolution  in  axial 
direction  and  in  torque  remain  unaffected  by  the  centrifugal  force 
field.  Thus,  equation  (25a)  becomes,  for  a  single  monotropic  mem¬ 
brane 


T  sin#  R  =  Tq  C  (4- A) 

Since,  in  general,  T/T ^  will  differ  from  unity,  the  fiber  paths  are 

obviously  no  longer  geodesics.  This  may  be  verified  by  comparing 
equation  (4- A)  with  Clairaut's  Theorem  given  in  the  form  of  equa¬ 
tion  (25a). 


The  axial  equilibrium  equation  (30)  becomes  for  a  single  mono¬ 
tropic  membrane: 

2  2  2 
n  T  cosjS  -  r  77  p  =  n  T  cosft  cosjS  -  R  r  77  p 


or 


V1  -  c2  -  !  (l-*2) 


=  cosa 


T  \2  C 


- P  (5-A) 

R 


where  K  is  defined  by 


K  = 


2 

277  r  p 
o 

n  T 


Solving  equation  (5-A)  for  cosOl  and  substituting  equation  (3-A)  for 
T/  Tq  yields  the  intrinsic  equation  for  the  meridional  shape  of  a 

monotropic  membrane  subject  to  both  a  uniform  internal  pressure 
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and  a  distributed  centrifugal  load  associated  with  a  rotation  around 
the  z-axis  and  a  uniform  mass  distribution  per  unit  length  of  the 
filaments : 

JT- 72  .  |  (i  -  r2) 

cosa  =  .  . . . . —  ( 6- A) 

Vii+“  (*-k2)]2  -  4 

R 

Thus,  the  equilibrium  shape  is  defined  by  three  characteristic 
parameters : 

(1)  The  parameter  defining  the  fiber  geometry:  C  . 

(2)  The  parameter  defining  the  pressure  loading:  K  . 

(3)  The  parameter  defining  the  centrifugal  loading:  £2  . 
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APPENDIX  B 


MERIDIONALLY  WOUND  ISOTENSOID  PRESSURE  VESSELS 

For  the  particular  case  of  C  -  0  (i.  e.  ,  meridional  fiber  di¬ 
rection  only)  and  for  p  =  p  =  Constant  ,  equations  (29)  and  (30) 
reduce  to 


r 


—  =  K  R 

PM 

(1-B) 

1  K,,  t,2v 

cosa  =  1  -  —  (1  -  R  ) 

(2-B) 

It  can  be  shown  (reference  18)  that  the  shape  of  a  buckled  slender 
column  with  large  deformations,  subject  to  eccentrical  loading, 
as  shown  in  figure  6,  satisfies  an  equation  of  the  form  (1-B)  as 
follows : 


Let  El  =  bending  stiffness  of  column 

-  P  =  compressive  load 

Then  the  equilibrium  condition  for  the  buckled  column  becomes 


^  El 

Pr  =  — 

P 


or 


El 


(3-  B) 


which  can  be  transformed  into  equation  (1-B)  by  setting 

P  r  2 

o 

K 


El 


Analytical  solutions  to  this  classical  problem  (MThe  Elastica") 
lead  to  tabulated  elliptical  integrals  of  the  first  and  second  kind  (ref¬ 
erence  13),  and  have  been  worked  out  in  detail  by  several  authors 
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(See,  for  instance,  references  19  and  20).  The  derivations  will  not 
be  reproduced  here,  but  only  the  principal  results  will  be  summa¬ 
rized. 


The  lengths  s  of  the  column  (or  meridional  lengths  of  the  iso- 

p  ro 

tensoid  pressure  vessel)  is  related  to  the  parameter  K  =  — — —  by 


(*,f) 


(4-B) 


where 


*  = 


3  :  Elliptical  integral  of  the  first  kind 

_  /k 


Characteristic  argument  of  elliptical  integral. 


The  maximum  deflection  r  of  the  column  (or  perimetral  radius  r 

o  o 

of  the  isotensoid  pressure  vessel)  is: 


r 

o 


(5-B) 


The  distance  h  between  ends  of  the  buckled  column  (or  between 
the  poles  of  the  isotensoid  membrane)  is: 


h_ 

r 

o 


(6-B) 


where 


£  :  Elliptical  integral  of  the  second  kind. 


The  "column  analogy11  for  the  meridional  isotensoid  uniform 
pressure  vessel  has  proven  useful  in  gaining  insight  into  the  types 
of  solutions  that  can  be  expected.  A  workable  analog  device  suitable 
for  general  study  has  been  constructed  from  a  thin  piano  wire,  loaded 
by  axial  forces  over  adjustable  bars,  as  shown  in  figure  7. 
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Several  solutions  for  the  case  C  =  0  are  shown  in  figure  11. 
Their  general  character  shows  a  peculiarity,  which  has  been  dis¬ 
cussed  in  reference  1. 

The  cusp  case  (K  =  2)  represents  a  closed  ovaloid  (Bermuda 
bottle).  Alternating  loops  crossing  the  z-axis  between  each  maxi- 
um  occur  between  K  =  2  and  K  =  4  .  These  loops  are  progres¬ 
sive  for  2  <  K  <  3.  32  and  regressive  for  3.  32  <  K  <  4  .  The 
transition  case  of  K  =  3.  32  represents  a  figure  resembling  a  lem- 
niscate  from  which  a  degenerated  closed  toroid  of  vanishing  inner 
perimeter  can  be  formed  by  rotation  around  the  z-axis. 

The  case  K  =  4  represents  an  aperiodic  loop  with  both 
branches  asymptotic  to  the  z-axis. 

For  K  >  4  ,  regular,  regressive  periodic  loops  such  as  shown 
in  figure  9(c)  occur. 
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APPENDIX  C 


SUMMARY  OF  DATA  REDUCTION  RELATIONS 


Calculation  of  Fiber  Tension  from  Pressure 

In  order  to  reduce  the  pressure- vs-elongation  data  to  terms  of 
tension  vs  strain  for  the  meridional  toroids,  the  effect  of  large  de¬ 
formations  must  be  taken  into  account.  This  is  accomplished  as 
follows : 

The  fiber  tension  in  the  meridional  winding  is  given  by: 


T  = 


2ir  p  r 

c 

n  K 
m 


where 


p  =  pressure 

r  =  radius  of  outer  equator 
o 

K  =  pressure  load  parameter 

n  =  number  of  fibers  in  the  meridional  windings 
m 

Let  denote  the  initial  value  of  the  outer  equatorial  radius,  at 

zero  pressure.  Let  e  denote  the  fiber  strain. 


Then, 


=  (On  (1  +  0 

o  \  of  o 


n  K 
m 


(ro)o2  <1+e> 
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Calculation  of  Volumetric  Strain 


The  volumetric  strain  is  defined  as  the  non-dimensional  linear 
scale  increase  which  would  correspond  to  a  given  increase  in  vol¬ 
ume  if  that  volume  increase  resulted  from  a  geometrically  similar 
expansion  of  the  pressure  vessel.  Thus,  if  is  the  "volumetric 

strain"  corresponding  to  a  volume  increase  AV  from  an  initial 

volume  V  ,  then, 
o 

vo  +  AV  =  Vo  (  +  V3 
This  yields  the  "volumetric  strain,  "  as, 


€ 


V 


-  1 
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TABLE  II 

DETAILS  OF  MEASUREMENTS 


MEASUREMENT 

INSTRUMENT 

PRECISION 

Meridian  length 

Extensometer 

0.  1% 

Outer  band  length 

Extensometer 

0.  1% 

Inside  diameter 

Steel  tape 

0.  1% 

Volume  change 

Platform  scale 

0.  05% 

(reservoir  wt.  ) 

Meridian  shape 

Polar- coord,  ring 
with  dial  gage  probe 

0.  1% 

of  min.  rad. 
of  curve 

Elev.  of  inner  equator 

Elev.  of  outer  equator 

► 

Height  gage  from 
horizontal  surface 

0.  2% 

Elev.  of  top  meridian 

to  marks  on  toroid 

Pressure 

Bourdon  gage  1 

0.  3% 
of  max. 

Meridians 

/ 


5 


Figure  5.  -  Pressure  Loaded  Shell  Made  from  Several 
Layers  of  Monotropic  Membranes. 
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Figure  8.  -  Isotensoid  Meridional  Shapes. 
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Figure  10.  -  Isotensoid  Domains  in  Function  of  the 
Parameters  C  and  K. 
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Figure  20.  -  Isotensoid  Shape  and  Filament  Geometry 
for  Closed  Torus  C  =  ±  .  35  Y.  =  4.  4. 
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Figure  24.  - 


Closed  Isotensoid  Pressure  Vessel  -  Type  IV. 


Figure  25.  -  Closed  Isotensoid  Pressure  Vessel  -  Type  V. 


Figure  26.  -  Typical  Cross  Section,  Meridionally- Wound  Toroid. 


Figure  27.  -  Toroid  Filament  Winder 


Polar  Coordinate  Measurement 


Figure  28.  -  Measurements  Made  on  Toroid. 


FIBER  TENSION  -  LBS 


3 


STRAIN  -  PERCENT 


Figure  30.  -  Typical  Load-Elongation  Data  For  Polyester  Fiber 


STRAIN  -  VS  -  TIME  MEASUREMENT 
SINGLE  FIBER 


Figure  31.  -  Strain-Time  Measurement  For  Single  Fiber 


Figure  33.  -  Two-L,obed  "Baseball  Seam"  Instability  of  Type  V 
Toroidal  Isotensoid  Pressure  Vessel. 


Figure  34.  -  Three-Lobed  "Baseball  Seam"  Instability  of  Type  V 
Toroidal  Isotensoid  Pressure  Vessel  (Top  View). 


Figure  35.  -  Thr ee-Lobed  "Baseball  Seam"  Instability  of  Type  V 
Toroidal  Isotensoid  Pressure  Vessel  (Side  View). 


Figure  36.  -  "Python"  Instability  of  Type  V  Toroidal  Isotensoid 
Pressure  Vessel. 


Figure  37.-  Fiber  Tension  vs.  Strain,  Toroid  #4 


FIBER  TENSION  -  LBS 
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STRAIN  -  PERCENT 


Center 

Top 

Bottom 


Figure  38.  -  Fiber  Tension  vs.  Strain,  Toroid  #5. 


Theoretical  Curve,  Initial  Size 
Measured  Data,  0.  5  Lb  Per  Fiber,  15  psi 
Theoretical  Curve  Expanded  By  7% 
Measured  Data,  9  Lb.  Per  Fiber  ,  240  ps.i 


Figure  39.  -  Cross-Section  Measurement,  Toroid  #5, 
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Theoretical  Curve  (From  Column  Analog) 
Measured  Data 


Figure  41.  -  Cross-Section  Measurement,  Toroid  #3. 
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mental  data,  obtained  from  filament -wound  toroidal  mental  data,  obtained  from  filament- wound  toroidal 

pressure  vessels,  confirm  the  validity  of  the  theory  pressure  vessels,  confirm  the  validity  of  the  theory 

and  demonstrate  application  of  the  analytical  design  and  demonstrate  application  of  the  analytical  design 

technique  for  filamentary  structures.  technique  for  filamentary  structures. 


